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Abstract 

Working at the Lifshitz level, we investigate the van der Waals interactions across 
a series of layers with a periodic motif. We derive the complete form of the van 
der Waals interaction as an explicit function of the number of periodic layers. We 
then compare our result with an approximation based on an anisotropic-continuum 
representation of the stratified medium. Satisfactory agreement between discrete- 
layer and continuum models is reached only for thicknesses of ten or more layers. 

1 Introduction 

Once one starts looking for them, multilayers are everywhere. (Bio) macro molecules such 
as phospholipid layers are prone to create multilamellar assemblies either in solution 
or in apposition to supporting interfaces [2]. Layer upon layer of polymer multilayers, 
assembled either electrostatically [3| or through interlayer hydrogen bonding [4 allow 
fabrication of multicomposite molecular assemblies of tailored architecture [5] with im- 
portant technological applications. Understanding molecular interactions in these systems 
is an important step in controlling the assembly process. Though these interactions are 
due to many different specific processes, the van der Waals interactions are their common 
feature. Unless one is satisfied with the trivial, pairwise additive formulation of van der 
Waals interactions in multilayer geometries, their general and exact derivation on the 
Lifshitz level is complicated, abstruse and seldom attempted. 
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Here we will set up a new approach to the van der Waals - Lifshitz interactions in 
multilayer geometries based on a recent reformulation of the Lifshitz theory jH] in terms 
of an algebra of 2x2 matrices that allows us to derive simple and transparent formulas for 
the van der Waals interactions across finely layered systems. We use this reformulation 
to evaluate the interactions across a multilayer and then compare these with a continuum 
uniaxial-dielectric approximation based on van der Waals interactions across anisotropic 
media. 



2 Formalism 

We consider a symmetric periodic array, Fig. ^ between a leftmost halfspace L and 
a rightmost halfspace R. The periodic motif is the sequence of AB pairs repeated N 
times between L and R, schematically LABAB...ABAR. Recently we showed [S] that in 
the Lifshitz theory computation of the secular determinant of the electromagnetic field 
modes can be mapped onto an algebra of 2 x 2 matrices. The secular determinant in 
fact follows from the 11 element of a transfer matrix that can be simply constructed from 
the interaction geometry as a product of discontinuity D and propagator T matrices. In 
the case considered here, following the "mnemonic" introduced in jH], the transfer matrix 
assumes the form 



x T A x B AB x T b x B ba x---xT a x B) ab xT b x B ba xT a x B al . (1) 



N 

In this notation the discontinuity and the propagator matrices become 

©ab = ( _^ ~ A ) = -Bba, (2) 

and 

T A,B = ( q e -2 PA , B a,b ^ , (3) 

where a and b are the thicknesses of the A and B regions, and 

A=( PA€B - pB€A ), (4) 

with ea{oj) and e B {uj) the frequency dependent dielectric functions of regions A and B. 
Finally 

P%b = Q 2 ~ £ ^f~. (5) 

c z 

We will assume that region A corresponds to hydrocarbon, £a{w), and regions L, B, R to 
water, e B (u), dielectric response. We use the standard [TU] forms for e^cj) and e B {u) 
where the dielectric response of water is described with one MW relaxation frequency, five 
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Figure 1: Symmetric periodic array between a leftmost region L and a rightmost region 
R. The periodic motif is the AB pair of layers repeated N times between L and R, 
symbolically LABAB...ABAR = L(AB) N AR. The thickness of layer A (hydrocarbon) is 
a, and that of layer B is b (water). Matrices D^, O^, D^^ and Ta,b ar e defined in the 
main text. 



IR relaxation frequencies and six UV relaxation frequencies and that of the hydrocarbons 
with four UV relaxation frequencies. 

The discontinuity matrix describes the propagation of the EM modes across the dielec- 
tric boundary and the propagator matrix their propagation inside a dielectrically homo- 
geneous region. The above equations are strictly valid for the TM field modes. The result 
appropriate for the TE field modes is obtained analogously, via a formal substitution 

\Pa£b + Pb£a) \PaPb + PbPaJ ' 

where pa,b are the magnetic permeabilities of the regions A, B. The transfer matrix can 
be written equivalently in the form 

M = B RA x A N xT A x B AL , (7) 
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Figure 2: Zero frequency n = van der Waals interaction free energy across a multilayered 
slab of layers A and layers B; a,b = 4nm (circles); a = lnm, b = 4nm (triangles); 
a = 8nm, b = 4nm (squares). Region A corresponds to lipid, regions B,L,R to water 
with static dielectric functions [TU] = 2 and €b,l,r = 80. The result of the exact Eq. 
I2H1 (bold curves) as well as approximate anisotropic continuum model Eq. EH (dashed 
curves) are presented. Clearly the exact and the approximate results differ only for small, 
N < 10, number of layers. 



where the matrix A can be obtained as 

T A x B AB xT B xD M 



1 - A. e~ 2pBb A (l - e - 2pBb ) 



The product A^ can be factored via the Abeles formula for square matrices [7j. This 
formula can be reproduced straightforwardly via induction starting from the trivial N = 2 
case |H] so that 



Q12 



sinh f \ sinh N£-fi2— sinh N£-?H sinh (N - l)f 

Vdet A V dct A 



(9) 



where 



^^(W 1 - 4 ^)- (10) 

Here and e - ^ are the two eigenvalues of the unitary matrix A* = A/ \J det A, viz. for 
unitary matrices the product of the eigenvalues equals one. From the definition of £, Eq. 
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Figure 3: Complete van der Waals interaction free energy Eq. 1^1 across N layers of 
hydrocarbon, A, and water, B, with thicknesses a,b = 4nm. The exact formula for N 
layers, Eq. EEl open circles, bold lines for the complete n sum, open squares, dashed line 
for the n = term (as on Fig. EJ). The approximate anisotropic continuum model, Eq. 
I3U1 with full circles, dashed lines for the complete n sum and full squares, dashed line for 
the n = term (same as on Fig. |2J). Exact and approximate results again differ only for 
small, N < 10, number of layers. 

ITTH we derive explicitly 

u 

£= (p A a + p B b) -log . = (p A a+pBb)-\ogf(u,p A a + p B b), (11) 

l + Jl- (ue-^+Psb)} 2 

where we defined 

2(1 "A') , n, I h\ U 

u— zzqz — and j(u,pAa+PBb) 



1 - A (e-M« + e~PB b ) + e ^( PA a+p B b) ' 1 + L _ ^ ue _ {pAa+PBb) ^ 

(12) 

The product A N can be further decomposed as 

fdet A) N / 2 pW-W 

= ' (i _ P. US) 

where 

F (N) = ( Pui ( Z\ P12( (at\ ) = ( 1 - e ~ 2N ") ~A= - e ~' C 1 - e-^" 1 ^ I, (14) 
V P2i(N) P22(N) J V ; v^letA v y 

with the unit matrix I. In Lifshitz theory the fluctuation free energy is directly related 
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Figure 4: The n = Hamaker coefficient, Eq. EI of the van der Waals interaction 
across a multilayered slab of hydrocarbon layers A and water layers B. Here a = lnm, 
b = 4nm (triangles), a,b = 4nm (circles) and a = 8nm, b = Anm (squares). Depending 
on the dielectric properties and the volume fractions of materials A and B the effect 
of the granularity persists for N on the order of 10 layers. Inset: the variation of the 
Hamaker coefficient at N = 5 layers as a function of the thickness of the lipid layer, a 
at a + b = 8nm. Clearly the largest Hamaker coefficient is obtained for larger volume 
fractions of the lipid material. 

to the 11 element of the tranfer matrix M [H]; here 

(detA) N / 2 e^ -1 ^ — — — 

m u = (1 _ e - 2 g) {pn(N)+p 21 (N)A RA + e- 2pAa A AL ( Pl2 (N)+p 22 (N)A RA )) . 

(15) 

The mode equation, or the secular determinant of the TM field modes can be written 

m™ = V TM {u, Q) = 0, (16) 

with an analogous equation for TE field modes. The combined secular determinant thus 
equals the product T)(u,Q) = T>tm(w,Q)Ute(w,Q). The free energy of the fluctuating 
EM modes in a system of iV (AB) layers can now be cast 6\ into a form containing the 
secular determinant of the TM and TE modes 

OO ' CO ' oo ' 

Q n=0 Q n=0 Q n=0 

(17) 

where u has been substituted by the imaginary Matsubara frequencies = i2n^n and 
the primed sum signifies that the n = term is taken with the weight 1/2. What we will 
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investigate now is the dependence of the van der Waals interaction free energy, defined as 
the difference 

F(iV; a, b) = F{N- a, b) - F{N — ► oo; a, b) (18) 

on the number of (AB) layers, N . From the above definition of the interaction free energy 
we obtain the corresponding Hamaker coefficient H(N; a, b) as 

3 Result 

We first write the 11 element for the TM modes in a form that makes explicit the depen- 
dence on N. We first notice that 



log mn = - log 2 sinh £+A log Vdet A e 5 +log (p n (N) + p 2 i(N)A RA + e' 2pAa A AL {p l2 {N) + p 22 (N)A RA )) . 

(20) 

so that the corresponding free energy Eq. El has three terms. The first does not depend 
on TV and is thus not of immediate concern here. The second term, linear in N, represents 
the chemical potential or the free energy of adding a single combined layer to the system. 
As already formulated in jH] this chemical potential has the form 

fj, = log VdetA A (+) , (21) 

where = is one of the eigenvalues of the matrix A*. However it is the last term of 
Eq. [201 that describes the interactions in the multilayer composite. Using Eqs. 1171 and IT81 
we end up with 

W ab)=kTTT log ( Pu(N) + P^ N ^ A + e- 2 ^ a A_AL(pn(N) +p 22 (AQA fM ; 

\Pii(oo)+p 21 (oo)A RA + e- 2 P^A AL (p 12 (oo)+p 22 (oc)A RA 

(22) 

the main result of this paper. It gives the complete van der Waals interaction across the 
stratified medium AB in a form with explicit dependence on the number of layers N. Not 
excessively complicated, it allows for straightforward numerical computations. 

Before analyzing Eq. E2in detail, we verify that it has the correct limiting behavior. 
When e A = e R , so that A = with p A = p R = p, if the stratified medium behaves 
homogeneously. In this case 

f(u, p A a + p B b) = 1 and £ = (p A a + p B b) = p(a + b). (23) 

Furthermore in the same limit we obtain for the matrix P 

P = 2 sinh {p A a + p B b) ( J e - 2N(p ° Aa+PBb) J • (24) 
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From these Eq. 12*2*1 becomes 

oo ' 

¥(N; a,b) = kTj2J2 lo § i 1 + ^RA^A L e- 2p{a+N{a+b)) ) . (25) 

Q n=0 

Thus we recover the standard Lifshitz expression for the interaction of media L and R 
across the medium A = B of thickness a + N(a + b). Obviously in the limit where the 
dielectric properties of A and B coincide, this is just the total thickness of the region A . 
The general formula Eq. E21 can be rewritten 

¥(N . aM = K f r QiQ l 0g ( gjW + MAQS M + e-^to.W + BB (AQS M ) 
27r ^ Jo \Pii(oo) +p 2 i(oo)A iM + e- 2 ^ a A AL (pi 2 (oo) +p 22 (oo)A iU ) 

(26) 

Its consequences can be fully appreciated only after a full numerical analysis. In Fig. 121 
we evaluate the n = term separately and then compare it to the full n summation that 
includes retardation effects, Fig. El 

4 Numerical computations and Conclusions 

For intuition about the general result Eq. |2Blwe evaluate the van der Waals interactions 
across a slab of material of thickness a+N(a+b), the continuum composite of layers A and 
B of Fig. ^ Because of its layered structure, we can associate with it a transverse dielectric 
function e± = e xx = e yy and a longitudinal dielectric function eii = e zz . The continuum 
composite is thus a uni-axial dielectric with two distinct values of the dielectric response 
parallel and perpendicular to the layer normal. In terms of the dielectric functions €a and 
6b by analogy to capacitors in series and in parallel, 

1 / \ , 1 1 f a b\ 

e± = —r (ae A + be B ) and - = —— — + — . 27) 
a + o e\\ a + o \6a £b J 

The longitudinal and the transverse dielectric response of a layered continuum composite 
thus depend on the volume fractions a/(a + b) and b/ (a + b) of the materials A and B in 
the system. 

We use the result of Sarlah and Zumer ^T] to derive the van der Waals interactions 
in a non-isotropic homogeneous uniaxial slab composed of periodic layers. If we associate 
the index A with the continuum uniaxial composite, we have for the TM modes 



The corresponding values for 



-V~ = \ Z A ~ + Te ) = ^ (29) 
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and mutatis mutandis ^T] for the TE field modes. Thus we obtain the interaction free 
energy as 

kT .°° 7 f°° 

¥(N; a, b) = / QdQ log (l + A m A^ L e~ 2 ^ a+N ^) , (30) 

n n=0 ^° 

plus an analogous term for the TE modes. This is the expression that we will compare 
with our exact formula Eq. [2H1 

Fig. El gives the dependence of the zero order term, n = 0, interaction free energy 
¥(N;a,b), in its exact form Eq. |2H]as well as its continuum approximate form Eq. 1201 
on the number of layers N, where we have taken A to correspond to hydrocarbon and 
B,L,R to correspond to water. The cases for b = 4nm and a = 1,4, 8nm are shown. It 
is clear from Fig. |2]that the effect of discreteness of layers is effectively gone after about 
10 layers in the stack. Also obviously the interaction free energy is larger the larger the 
volume fraction of the lipid in the system. 

The same is true also for the complete form, Eq. EBl with the explicit summation over 
n, Fig. El The higher order terms in the n-summation are obviously important only for 
small N. In both cases, n = as well as full summation over n, the way the continuum 
limit is approached obviously depends on the characteristics of the layers, especially on 
their respective thicknesses. As is well known ^U] for the hydrocarbon-water system the 
n = term gives the dominant contribution to the interaction energy and assures that 
contrary to the case of other materials the van der Waals interactionshave a non-retarded 
form even at very large separations. 

The n = Hamaker coefficient, Fig. 0] and Eq. EH reinforces conclusions based on 
the evaluation of the interaction free energy. Again the continuum limit of the Hamaker 
coefficient is safely reached after about 10 layers, the swiftness of the transition to the 
continuum approximation being dictated by the thickness and the dielectric properties of 
the layers A and B. Also the magnitude of the Hamaker coefficient, for fixed thickness of 
the periodic layer, a + b = 8nm in the case examined, depends on the volume fraction of 
both layers. It is biggest when the lipid layer volume fraction is highest. This is indeed 
as one would expect since for high water volume fractions the dielectric inhomogeneity 
and thus van der Waals interactions are non-existent. It is only at much higher values 
of N that the effect of retardation will set in, leading to a levelling off of the Hamaker 
coefficient at about half of its n = value. Since our focus here is the granularity or 
stratification of the medium, we do not delve any further on the retardation effects. 

If the Hamaker coefficient were evaluated in the model based on an isotropic mixture 
of A and B the Hamaker coefficient would be about 20% smaller than in the anisotropic 
continuum model. 

Comparing the exact dependence of the van der Waals interactions across a periodic 
multilayer on the number of layers with a continuum approximation based on the van der 
Waals interactions across uni-axial anisotropic composite, we have been able to assess the 
validity of the continuum approximation and show that it indeed reduces to the properly 
defined limit after about 10 layers. Again, this limit has to acknowledge the uniaxial 
continuum structure of the multilayer. 
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